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ABSTRACT 

We calculate neutrino cross sections from neutral current reactions in the 
dense matter encountered in the evolution of a newly born neutron star. Effects 
of composition and of strong interactions in the deleptonization and cooling 
phases of the evolution are studied. The influence of the possible presence of 
strangeness-rich hyperons on the neutrino scattering cross sections is explored. 
Due to the large vector couplings of the E~ and S~, \Cv\ ~ 2, these particles, if 
present in protoneutron star matter, give significant contributions to neutrino 
scattering. In the deleptonization phase, the presence of strangeness leads to 
large neutrino energies, which results in large enhancements in the cross sections 
compared to those in matter with nucleons only. In the cooling phase, in which 
matter is nearly neutrino free, the response of the Tr hyperons to thermal 
neutrinos is the most significant. Neutrinos couple relatively weakly to the A 
hyperons and, hence, their contributions are significant only at high density. 

Subject headings: dense matter - stars: neutron - stars: opacities - stars: 
neutrinos 
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1. INTRODUCTION 

The general nature of the neutrino signature expected from a newly formed neutron 
star (hereafter referred to as a protoneutron star) has been theoretically predicted (Burrows 
& Lattimer 1986) and confirmed by the observations (Bionta et al. 1987; Hirata et al 1987) 
from supernova SN1987A. Although neutrinos interact weakly with matter, the high baryon 
densities and neutrino energies achieved after the gravitational collapse of a massive star 
(> 8 solar masses) cause the neutrinos to become trapped on the dynamical time scales 
of collapse (Sato 1975; Mazurek 1975). Trapped neutrinos at the star's core have Fermi 
energies E^, ~ 200 — 300 MeV and are primarily of the Ue type. Thermally produced 
and r neutrino pairs are also trapped, but with zero chemical potential. Neutrinos escape 
after diffusing through the star exchanging energy with the ambient matter, which has an 
entropy per baryon of order unity in units of Boltzmann's constant. Eventually they emerge 
from the star with an average energy ~ 10 — 20 MeV and in nearly equal abundances of all 
three flavors, both particle and anti-particle. 

Neutrino interactions in dense matter have been investigated by various authors 
(Tubbs & Schramm 1975; Sawyer 1975,89,95; Lamb & Pethick 1976; Lamb 1978; Sawyer 
& Soni 1979; Iwamoto & Pethick 1982; Iwamoto 1982; Goodwin &Pethick 1982; Burrows 
& Mazurek 1982; Bruenn 1985; van den Horn & Cooperstein 1986; Cooperstein 1988; 
Burrows 1988; Horowitz and Wehrberger 1991a,b;1992; Reddy & Prakash 1995). The 
charged current absorption and neutral current scattering reactions are both important 
sources of opacity. The neutral current scattering involves all flavors of neutrinos scattering 
on nucleons and leptons. Scattering from electrons is important for energy and momentum 
transfer (Tubbs and Schramm, 1975). However, for lepton number transport, nucleon 
scattering and absorption are the dominant processes. 

Surprisingly little attention has been paid to the effects of composition and of strong 
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interactions of the ambient matter on neutrino opacities. In the few attempts to date, the 
effect of interactions was investigated for non-degenerate nuclear matter by Sawyer (1975; 
1989) and for degenerate pure neutron matter by Iwamoto & Pethick (1982). Treating 
nucleons in the non-relativistic hmit, these calculations predict an increase in the mean free 
path by a factor of ~ 2 — 3, for (2-4) times the nuclear density. More recently, relativistic 
calculations based on effective Lagrangian models for dense neutron star matter have been 
performed by Horowitz and Wehrberger (1991a,b;1992). Here, the differential cross-sections 
for matter containing nucleons and electrons were calculated using linear response theory. 
A reduction of 30-50% over the case of non-interacting nucleons was reported in these 
calculations. In addition to strong interaction modifications, electromagnetic interactions 
can increase the mean free path by 50-60% for electron type neutrinos through collective 
effects (Horowitz 1992). So far, the influence of interactions has been investigated in 
protoneutron star calculations only by a simple scaling of the non-interacting results 
(Burrows 1990, Keil 1994). Furthermore, there have been no calculations performed 
including the multi-component nature of the system. We note that Keil and Janka (1995) 
have recently carried out deleptonization and cooling simulations including hyperons in the 
EOS, but they ignored opacity modifications. We view it as essential that opacities be 
consistent with the composition, which has not been a feature of protoneutron star models 
to date. The composition of the star changes significantly from the deleptonization phase, 
in which neutrinos are trapped, to the cooling phase, in which only thermally produced 
neutrinos are present. So far, important opacity modifications due to the changing lepton 
content and composition have not been treated satisfactorily in evolutionary calculations. 

Although the composition and the equation of state (EOS) of the hot protoneutron 
star matter are not yet known with certainty, QCD based effective Lagrangians have opened 
up intriguing possibilities (Kaplan & Nelson (1986); Glendenning 1986,1992; Glendenning 
& Moszkowski (1991); Kapusta & Ohve (1990); Ellis, Knorren & Prakash (1995); Knorren, 



- 5 - 



Prakash & Ellis (1995), Prakash, Cooke & Lattimer 1995). Among these is the possible 
existence of matter with strangeness to baryon ratio of order unity. Strangeness may be 
present either in the form of fermions, notably the A and E~ hyperons, or, in the form 
of a Bose condensate, such as a K"- meson condensate, or, in the form of s quarks. In 
the absence of trapped neutrinos, strange particles are expected to appear around 2 — 4 
times the nuclear matter density of Uq = 0.16 fm~^. Neutrino-trapping causes the strange 
particles to appear at somewhat higher densities, since the relevant chemical potential 
II — lie — IJ'Ve matter with high lepton content is much smaller than in the untrapped case 
(Thorsson, Prakash & Lattimer 1994; Ellis, Knorren & Prakash 1995; Knorren, Prakash & 
Elhs 1995, Prakash et al 1996). 

A new feature that we consider here is the role of strangeness. To date, only neutrino 
opacities for strange quark matter have been calculated (Iwamoto, 1982). Here, we study 
neutrino scattering mean free paths in matter containing strangeness in the form of 
hyperons. Specifically, we calculate neutrino opacities from neutral current reactions in 
matter containing hyperons and which are faithful to the EOS. In a first effort, this will 
be achieved using a mean field theoretical description which includes hyperonic degrees of 
freedom. This approach has several merits. For example, aspects of relativity, which may 
become important at high density, are naturally incorporated. Modifications of the opacity 
due to correlations (RPA) are also possible in such an approach. Further, comparisons 
with alternative potential model approaches (Iwamoto & Pethick, 1982; Sawyer, 1989) are 
straightforward. Neutrino opacities in matter containing other forms of strangeness will be 
considered in a separate work. Contributions from charged current reactions are essential for 
a complete description of the protoneutron star evolution. With appropriate modifications 
of the formalism presented in this work, calculations that include the compositional changes 
in the distinct phases of the evolution will be reported in a later work. 
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In §2, the formalism to calculate the neutrino scattering cross sections in a multi- 
component system is discussed. New analytical formulae for the response functions, which 
facilitate accurate calculations of the cross sections in all regimes of matter degeneracy 
are derived. In §3, the composition of beta-equilibrated matter with and without strange 
baryons is determined based on a field theoretical model. This section also contains a 
description of the relevant physical conditions in the evolution of a newly born neutron 
star. In particular, the composition in the distinct phases of deleptonization and cooling 
are discussed. Results of the neutrino scattering cross sections in these two phases are 
presented in §4. Conclusions and directions for future study are given in §5. 



2. NEUTRINO INTERACTIONS WITH BARYONS 

Neutrino interactions with matter proceed via charged and neutral current reactions. 
The neutral current processes contribute to elastic scattering, and charged current reactions 
result in neutrino absorption. The formalism to calculate neutral current scattering rates 
in dense matter is summarized below. The interaction Lagrangian for neutrino scattering 
reactions is based on the Wienberg-Salam-Glashow theory (Weinberg 1967; Salam 1968; 
Glashow 1961): 

= for ^ + B ^ + B , (1) 

where Gp ^ 1.436 x 10^"^^ erg cm~^ is the weak coupling constant. The neutrino and target 
particle weak neutral currents appearing above are: 

= ^,r(Cvi-CAa,)^l^i, (2) 

where i — A, E~, E+, E°, S~, • • • and e~,ii~. The neutral current process couples 
neutrinos of all types (e, /i, and r) to the weak neutral hadronic current j^. The vector 
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and axial vector coupling constants, Cvi and Cai are listed in Table 1. Numerical values 
of the parameters that best fit data on charged current semi-leptonic decays of hyperons 
are (Gaillard & Sauvage 1984): D=0.756 , F=0.477, sin2^vK=0.23 and sin^c = 0.231. Due 
to the large vector couphngs of the T,~~ and S~, |Cy| ~ 2, these particles, if present in 
protoneutron star matter, will give significant contributions to neutrino scattering. Neutrino 
scattering off leptons in the same family involves charged current couplings as well, and one 
has to sum over both the contributing diagrams. At tree level, however, one can express the 
total coupling by means of a Fierz transformation; this is accounted for in Table 1. 

Given the general structure of the neutrino coupling to matter, the differential cross 
section for elastic scattering for incoming neutrino energy E^, and outgoing neutrino energy 
El is given by (Fetter and Walecka 1971) 

where the neutrino tensor La/s and the target particle polarization 11"'^ are 

Lap = S[2kakp + {k ■ q)gai3 - {kaQp + qakp) T ieaptivk^q''] , (4) 
= -iJ^,TT[G\p)JaG\p + q)M. (5) 

Above, k/^ is the incoming neutrino four momentum and is the four momentum transfer. 
The Greens' functions G^{p) (the index i labels particle species) depend on the Fermi 
momentum kpi of target particles. In the Hartrcc approximation, the propagators are 
obtained by replacing Mj and kpi in the free particle propagators by M* and kp^ (see 
below), respectively. The current operator is 7^ for the vector current and 7^75 for the 
axial current. Given the V-A structure of the particle currents, we have 

For the vector polarization, { Ja, Jg} :: {ja^lf}}, ior the axial polarization, 
{Ja,Jf3} ■■ {7a75, 7/375} and for the mixed part, {Ja, Jp} " {7a75,7^}- Further, the 
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polarizations contain two functions, the density dependent part that describes particle-hole 
excitations and the Feynman part that describes particle-antiparticle excitations. For elastic 
scattering, with < 0, the contribution of the Feynman parts vanish. Using vector current 
conservation and translational invariance, may be written in terms of two independent 
components. In a frame where = {qo, \q\, 0, 0), we have 



= n^2 and = ■ 



The axial current-current correlation function can be written as a vector piece plus a 
correction term: 

K = K + 9,uli^ ■ (7) 
The mixed, axial current-vector current correlation function is 

= ie^,.,a,o?"n^^ . (8) 

The above mean field or Hartree polarizations, which characterize the medium response 
to the neutrino, have been explicitly evaluated in previous work (Horowitz & Wehrberger 
1991). In terms of these polarizations, the differential cross-section is 

with 

2A;o(A;o - gp) + 

A = —r^ ^ — ; B = 2ko-qo. (10) 

The polarizations may be combined into three uncorrelated response functions i?i,i?2 and 
i?3 by summing over the contributions from each particle species i: 

Ri = J:[Cv^ + Cl][ImUl, + lmW^], (11) 

i 

R2 ^ Cvi Im + Cl [Im - Im , (12) 

i 

i?3 = ±E 2C^,C^,ImnV^. (13) 
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The imaginary parts of the lowest order one-loop polarization parts, required for the 
differential cross-sections, are evaluated at finite density and temperature. These retarded 
polarizations characterize the medium response at the mean field level. These functions 
depend upon the individual concentrations, which are controlled by the effective chemical 
potentials Vi and the corresponding effective masses M*, for which a many-body description 
of the multi-component system is required (see §3). 

Analytic closed form expressions for the density dependent polarization functions at 
zero temperature have been evaluated by Lim and Horowitz (1989). The finite temperature 
polarizations have been investigated by Saito et al. (1989). However, we are not aware that 
closed form analytical formulae for the finite temperature retarded polarizations have been 
been reported elsewhere in the literature. Here, we provide simple analytical formulae for 
these polarizations in the space like region (go < 1^1 and < 0). The integral forms of the 
finite temperature correlation functions and their simplification are discussed in Appendix 
A. The results for the various polarizations (for one species of the target particles) are: 

2 ■ 4g2 



2 

ImHr = -\[/2 + go/i]+ ^ 



Im Hi 



(14) 
(15) 



ImH^ = -—lo, (16) 
Znq 

ImUvA = ^[go/o + 2/i]. (17) 

Above, Io,h, and I2 are phase space integrals. Each phase space integral consists of two 
contributions: 

+ i = 0,l and 2, (18) 

where the superscript + refers to particle excitations and the superscript — to antiparticle 
excitations. At finite temperature, both particle and antiparticle excitations contribute, 
although the latter contribution is exponentially suppressed by the factor exp(— ///T). 
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Explicit analytical forms for the phase space integrals are: 

= qo-nt: (19) 
It = ±uqo-^-T'Ct-e.T^t, (20) 

+ 2T^ii-2e.T^Ct + e'_Tit, (21) 
where the factors ^t be expressed as differences of polylogarithmic functions Lin as 

Ct = Lzr,{-af) - Lzn{-ai) ■ (22) 



= exp(^ — j and a2 = exp j . (23) 

The polylogarithimic functions 

Li„(z) = r :^!!iziM (i^; , Lii{x) = log(l - x) (24) 
are defined to conform to the definitions of Lewin (1983). 

The factor e_ in the above equations is the kinematic lower limit and is given by 



qo , q 



M*2 

v^v- '''' 



The quantities M* and v are the effective mass and effective chemical potential of the 
particle, respectively (see §3). 

The above formulae are vahd for go > 0. We can extend them to the case qo < using 
the principle of detailed balance, which ensures that 

Im U{qo < 0) = exp(-|go|/T) Im U{qo > 0) . (26) 



These analytical expressions allow us to understand the behavior of the correlation 
functions with both density and temperature. This particular representation for the 
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relativistic response functions permits one to calculate accurately the cross sections in 
all regimes of matter degeneracy. This is especially important in multi-component dense 
matter, where one expects to encounter different levels of degeneracy for each particle 
species. Note also that this approach incorporates the effects of strong interactions through 
modifications of the particle propagators. In addition, the time consuming and often 
problematic numerical integration required to evaluate the phase space integrals (see 
Eq. ( P^ in Appendix A) is avoided. 

The total inclusive scattering rate for neutrinos from a hot and dense system is 

obtained by integrating over the allowed kinematic region in the uj — \q\ space. Explicitly, 

a{E) r2E-qo 
V = 2^ i-oo L\ ^"^^^ ^^"^ ^^'^ ■ ^^^^ 

Note that the total cross-section per unit volume, which has the dimension of inverse length, 

gives the inverse collision mean free path. Various other transport coefficients, such as the 

diffusion coefficient and the thermal conductivity, can also be obtained in a similar fashion. 



3. COMPOSITION OF NEUTRON STAR MATTER 

To explore the influence of the presence of hyperons in dense matter, we employ a 
relativistic field theoretical model in which the interactions between baryons are mediated 
by the exchange of (t,uj, and p mesons. The full Lagrangian density is given by (Serot & 
Walecka 1986), 

L = Lh + Li 

B 
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I 

Here, B are the Dirac spinors for baryons and t is the isospin operator. The sums include 
baryons B = n,p, A, S, and S, and leptons, / = e~ and fi~ . The field strength tensors for 
the u and p mesons are Wfj,i, = dfj,ujy — d^uj^ and 'Q^u = d^h^, — di^h^ , respectively. The 
potential U (a) represents the self-interactions of the scalar field and is taken to be of the 
form 

Uia) = ^bM^ig^r^af + ^cig^N^)^ ■ (28) 

Electrons and muons are included in the model as non-interacting particles, since their 
interactions give small contributions compared to those of their free Fermi gas parts. 

In the mean field approximation, the partition function (denoted by Zh) for the 
hadronic degrees of freedom is given by 

In Zh = (3V [»o' + >?&o - Wy " ^(^) 

where (3 = {kT)~^ and V is the volume. The contribution of antibaryons is not significant 
for the thermodynamics of interest here, and is therefore not included in Eq. (E^). Here, the 



effective baryon masses = Mb — gasO' and E*^ = ^k? + Mff. The chemical potentials 
are given by 

PB = J^B + QloB^O + gpBtsBbo , (30) 

where t^s is the third component of isospin for the baryon. Note that particles with = 0, 
such as the A and T,^ do not couple to the p. The effective chemical potential ub sets the 
scale of the temperature dependence of the thermodynamical functions. 

Using Zh, the thermodynamic quantities can be obtained in the standard way. The 
pressure Ph = TV~^ ^^Zh, the number density for species B, and the energy density Eh 
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are given by 



e„ = \mW + U(o) + \ml^l + im>S + 2 E / (^i^J (eW^i-' + l)'' , (31) 
The entropy density is then given by sh = P{£h + Ph — J2b f^B^B)- 

The meson fields are obtained by extremization of the partition function, which yields 
the equations 



m: 



B B 

73 



^^E«..^/^f(e— .1)-'. (3.) 



B 

The total partition function Ztotai = ZhZl, where Zl is the standard non-interacting 
partition function of the leptons. 

The additional conditions needed to obtain a solution are provided by the charge 
neutrality requirement, and, when neutrinos are not trapped, the set of equilibrium chemical 
potential relations required by the general condition 

Aij = bif^n - Qiii'i , (33) 

where hi is the baryon number of particle i and qi is its charge. For example, when i = e~ , 
this implies the equalities 

/iA = ^J'T,o = fiBO = /in , 
A^E- = /^H- = /in + /ie , 

/ip = /iS+ = fJ'n- f^e- (34) 

In the case that the neutrinos are trapped, Eq. (|33|) is replaced by 

IJ'i = bijJn - qiif^l - jJ'uJ ■ (35) 



-14- 



The new equalities are then obtained by the replacement /ig — /ie — fiue iii Eq. (pi]). The 
introduction of additional variables, the neutrino chemical potentials, requires additional 
constraints, which we supply by fixing the lepton fractions, Yie, appropriate for conditions 
prevailing in the evolution of the protoneutron star. The contribution to pressure from 
neutrinos of a given species is P^, = (l/247r^)/i^. 

In the nucleon sector, the constants g^N, g^N, gpN,b, and c are determined by 
reproducing the nuclear matter equilibrium density Uq = 0.16 fm~^, and the binding energy 
per nucleon (~ 16 MeV), the symmetry energy (~ 30 — 35 MeV), the compression modulus 
(200 MeV <Ko< 300 MeV), and the nucleon Dirac effective mass M* = (0.6 - 0.7) x 939 
MeV at riQ. Numerical values of the coupling constants so chosen are shown in Table 2. 
This particular choice of model parameters are from Glendenning and Moszkowski (1991) 
and will be referred to as GMl hereafter. The prevalent uncertainty in the nuclear matter 
compression modulus and the effective mass M* does not allow for a unique choice of these 
coupling constants. The high density behaviour of the EOS is sensitive to the strength of 
the meson coupling constants employed. Lacking definitive experimental and theoretical 
constraints, this choice of parameters may be considered typical. 

The hyperon coupling constants may be determined by reproducing the binding energy 
of the A hyperon in nuclear matter (Glendenning & Moszkowski 1991). Parametrizing the 
hyperon-meson couplings in terms of nucleon-meson couplings through 

= gan/gaN, x^jH = goju/guiN, = gpu/ gpN 1 (36) 
the A binding energy at nuclear density is given by 

(5/A)a = -28 = x^A^f^ATtUo - x^^g^NCTo , (37) 

in units of MeV. Thus, a particular choice of x^a determines x^^a uniquely. To keep the 
number of parameters small, the coupling constant ratios for all the different hyperons are 
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assumed to be the same. That is 



*^(T "^(tK "^CrS •^CF^ 0.6 , 



(38) 



and similarly for the a; 




0.653. 



(39) 



The p-coupling is of less consequence and is taken to be of similar order, i.e. Xp = . 



the relative fractions, Yi = Ui/nb, of the baryons and leptons in charge neutral and 
/3— equilibrated neutrino-free matter at zero temperature. The upper panel shows the 
concentrations for the case in which only the nucleonic degrees of freedom are allowed. The 
lower panel contains results for the case in which hyperonic degrees of freedom are also 
allowed. For the model parameters chosen, the S~ hyperon appears at a density lower than 
the A hyperon. This is because the somewhat higher mass of the S~ is compensated by the 
presence of the e~ chemical potential in the equilibrium condition of the S~. More massive 
and more positively charged particles appear at higher densities. With the appearance of 
the negatively charged which competes with leptons in maintaining charge neutrality, 
the lepton concentrations begin to fall. The important point is that, with increasing density, 
the system contains many baryon species with nearly equal concentrations. 

It should be pointed out, however, that moderate changes in the poorly known E and 
H couplings has large effects on the appearance of negatively charged particles (Prakash et 
al., 1996). Increasing the coupling constants of a hyperon species delays its appearance to 
a higher density. This is because the threshold condition, Eq. (pSf ), receives contributions 



3.1. Composition in a Cold Catalyzed Neutron Star 



Old and cold neutron stars are essentially neutrino free. In Figure we show 
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from the cr, a;, and p fields, the net result being positive due to that of the uj. If all the 
couplings are scaled up, the positive contribution becomes larger, and hence the appearance 
of the particle is delayed to a higher density. Clearly, the thresholds for the strange particles 
are sensitive to the coupling constants, which are presently poorly constrained by either 
theory or experiment. Notwithstanding these caveats, it is clear that one or the other 
hyperon species is likely to exist in dense matter. 

3.2. Composition in an Evolving Newly Born Neutron Star 

The protoneutron star formed subsequent to the core bounce, and its early evolution 
has been investigated in earlier works (Burrows & Lattimer 1986, Keil & Janka 1995). 
Detailed studies of the dynamics of core collapse and supernova indicate that within 
milliseconds of the shock wave formation the core settles into hydrostatic equilibrium, with 
a low entropy and large lepton content. The electron lepton fraction Yl^ — Ye + at 
bounce in the interior is estimated to be about 0.4. The electron type neutrinos formed 
and trapped in the core during collapse are degenerate with a chemical potential of about 
300 MeV. In the case of muons it is generally true that, unless fi > m^c^, the net number 
of /x's or i/^'s present is zero. Because no muon-flavor leptons are present at the onset of 
trapping, Y^^ = —Y^. Following deleptonization, Y^^ = 0, and is determined by //^ = /le 
for lie > m^c^ and is zero otherwise. 

The entropy per baryon, in the interior is low, of order unity, which corresponds to 
temperatures between 5-30 MeV in the interior. The temperature and entropy increases as 
one moves from the center outwards (Burrows & Lattimer 1986). At the very early stage, 
the electron neutrino chemical potential is largest in the center and drops appreciably as a 
function of the distance from the center. This gradient in the electron neutrino chemical 
potential is primarily responsible for driving the deleptonization phase. During the early 
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deleptonization phase, the neutrino re-heating of matter increases the central temperature. 
This, coupled with other diffusive processes, reverses the temperature gradient. On time 
scales of about 10-15 s, the central temperature is raised to about 50 MeV {S ~ 2) and 
decreases outwards. This marks the onset of the cooling phase in the central regions of 
the star. The conditions at the onset of deleptonization and cooling are thus significantly 
different. In summary, the composition and temperature in the central regions of the star 
at the beginning of deleptonization are characterized by a high lepton fraction {Yie = 0.4) 
and low entropy {S ~ 1), while the cooling phase is characterized by a low neutrino fraction 
(y^/g ~ 0) and high entropy (5 ~ 2). 

For a full treatment of neutrino transport during the evolution, opacities for a wide 
range of composition and matter degeneracy are required. The dynamical changes in the 
lepton fraction and the temperature modify the composition of matter and the typical 
neutrino energies in the inner core. In addition, there are structural changes in the interior 
associated with changes in Yie and T (Prakash et al., 1996). Computer simulations of the 
evolution account for these effects dynamically. Here, we have identified two important and 
distinct phases during the early evolution based on such simulations (Burrows & Lattimer 
1986; Keil & Janka 1995), in order to highlight the role of neutrino scattering. 

3.2.1. The Deleptonization Phase 

For typical conditions in this phase, Yl^ = 0.4 and T = 20 MeV, Figure |^ shows 
the concentrations of the various species. The top panel shows results for matter with 
nucleons only. The bottom panel refers to the case in which hyperons are present. In both 
cases, neutrino trapping significantly alters the composition from the neutrino-free case 
(see Figure |lD. This is because fi = fie — ^iv^ is much smaller than /ig in the neutrino- free 
case. This results in large electron, and hence, to satisfy charge neutrality, large proton 
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concentrations. Also, the appearance of hyperonic components is delayed to higher 
densities; in particular, the concentrations of electrically charged hyperons are suppressed. 
Finite temperature effects on the composition are less significant, and in general favors 
the presence of strange baryons. Relative to nucleons-only matter, the neutrino chemical 
potential and hence the neutrino concentrations increase substantially with density in 
matter containing hyperons. This introduces an important distinction between the opacities 
in matter containing only nucleons and in matter containing hyperons as well. 

3.2.2. The Coolmg Phase 

In this phase, ~ 0. Figure |^ shows the concentrations Yi for a typical T = 50 MeV. 
Here, the strangeness-bearing components occur at significantly lower densities than in 
the neutrino-trapped case. A comparison of the two panels shows that in strangeness-rich 
matter, the S~ hyperons effectively replace the leptons in maintaining charge neutrality. 
The increasing abundance of the neutral particles A and TP has important consequences for 
neutrino scattering. Although the lowest order contributions from the neutral A and TP are 
smaller than those of the T^ and (see Table 1), their presence furnishes baryon number, 
which decreases the relative concentrations of nucleons. (In contrast, neutrons are the most 
abundant particles in nucleons-only matter.) Relative to neutrino trapped matter (Figure 
the larger temperature and the absence of a neutrino chemical potential both contribute 
synergetically to enhance the hyperonic fractions. 

There exists an unambiguous difference between the deleptonization (neutrino-trapped) 
and cooling (neutrino-free) phases (Prakash et al., 1996). Neutrino trapping delays 
the appearance of hyperonic components to higher baryon densities. This implies that 
during the early deleptonization phase, matter consists mostly of non-strange baryons, 
except possibly at high densities. The cooling phase is characterized by the presence of a 
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substantial amount of strangeness-rich hyperons, as they appear at lower densities. Thus, 
we may expect modifications to the neutral current scattering due to strangeness to be 
quantitatively different during the two phases. In the following section, we present the 
scattering cross sections in both these phases. 

4. NEUTRINO SCATTERING CROSS SECTIONS 

We turn now to quantitative results for the neutrino scattering rates. The differential 
and total cross sections are evaluated per unit volume of matter, and the contribution from 
each particle species is summed over. The results presented here take into account strong 
interaction and finite temperature effects through the finite temperature mean field response 
functions. In order to calculate the differential cross sections using Eq. (P), the chemical 
potentials of all particle species, the temperature and the mean fields need to be specified. 
These are provided by the calculations described in §3. For the EOS employed, the density 
in the central regions of the star is in the range of 3 — 7 times the nuclear saturation density 
(Prakash et al 1996), so for the most part, we choose to focus on a representative density of 
ub = 0.64 fm"^. 

4.1. Influence of Composition and Temperature 

A first orientation to the interaction corrections and finite temperature effects on the 
neutrino scattering rates is provided by investigating pure neutron matter. Figure ^ shows 
the effect of varying density on the differential cross-sections (solid curves 1(a) and 2(a)). 
The curve labels 1 and 2 refer to ub = 0.32 fm~^ and ub = 0.64 fm"'^, respectively. For 
comparison, the corresponding results for a free Fermi gas are also shown (dashed curves 
1(b) and 2(b)). Relativistic effects from the matrix elements and phase space considerations 
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are both small for non-interacting baryons. These effects, however, significantly influence 
the mean field response. In particular, the response is sensitive to the effective mass and the 
effective chemical potential. Scalar interactions, which generate nucleon effective masses 
that decrease with increasing density, considerably alter the response in comparison to that 
of a free Fermi gas. For small u;, interactions lead to a suppressed cross section, while the 
large uj cross sections are enhanced. 

Figure |^ illustrates finite temperature effects on the differential cross sections in 
pure neutron matter. At zero temperature, only the positive uj excitations are present, 
since negative uj excitations are blocked by the Pauli principle. However, changing the 
temperature on the scale of the energy transfer uj furnishes sufficient number of particles in 
excited states. This gives rise to a significant amount of negative uj response (Iwamoto & 
Pethick 1982; Sawyer 1989). This is true even if the temperature is very small compared 
with the target particle Fermi energies. For positive uj, finite temperature corrections 
depend on the density. For low baryon densities, even temperatures of order a few tens of 
MeV render the system non-degenerate. Hence, the positive uj response shows considerable 
sensitivity to temperature (upper panel). With increasing density, when matter becomes 
increasingly degenerate, effects of temperature are negligibly small for > (lower panel). 

The response of a multi-component system is significantly different from that of a 
single component system. At a given density, the concentrations and the effective masses 
of the individual particles, and, their specific coupling to the neutrinos determine the total 
response. These features are illustrated in Figures ^ and |^, where we contrast the differential 
cross sections aX ub = 0.64 fm~^ and T = 30 MeV in nucleonic matter and in matter with 
strange baryons. The contribution of each particle species is shown separately in these 
figures. At zero temperature (Figure |^), the particles are all degenerate and only positive 
UJ excitations are present. The contributions from the S~ and the A dominate the neutron 
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contribution in the low uo region (see lower left panel). This may be roughly understood 
in terms of the magnitudes of the product of two factors for the various particles: The 
first factor is = + Mf (true only at T = 0), which reflects the composition and the 
effective mass. The second factor is Cy_ + 3C^. and is a measure of the neutrino coupling. 
The dominance of the S~ is due to the fact that both factors above are much larger 
than those for the A and the neutron. Notice also that since z/a > z/„ at this density, the 
contribution from the As is comparable to those from neutrons despite the fact its effective 
coupling is relatively smaller (see Table. 1). At finite temperature (Figure the uj < Q 
response is opened up. Further, the sharp fall off in the cross section (for > 0), associated 
with the S~ hyperon, is smoothed out. Due to the varying degeneracies, the response of the 
different particle species are different. Most significant changes occur in the region where 
uj/T ~ 1. 

In what follows, we study the differential and total cross sections at finite temperatures, 
densities and composition relevant to the early evolution of a newly born neutron star. 

4.2. The Deleptonization Phase 

Here, we calculate the cross sections at a baryon number density ub = 0.64 fm~^, 
temperature T = 20 MeV and lepton fraction Y^g = 0.4, which are representative values in 
the central regions at the beginning of this phase. The composition under these conditions 
was discussed in §3.2.1. Since neutrinos are degenerate, the relevant neutrino energies lie in 
the range 250-450 MeV. The effects of neutrino degeneracy are incorporated by explicitly 
by including the neutrino final state blocking factor 1 — f{E — fi^,^ + u) in the expression 
for the differential scattering cross section. The neutrino chemical potential in the interior 
increases with density (since Yie is held constant at the value 0.4). At a given density, only 
neutrinos close to the Fermi surface can actively participate in the diffusion process. The 
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different neutrino energies chosen reflect the different neutrino chemical potentials in matter 
with nucleons and hyperonic components. Figure ^ shows a comparison of the differential 
cross sections in normal and strangeness-rich matter. Note that electrons, due to their large 
concentrations, contribute nearly as much to the cross sections as the neutral particles. 
Further, the hyperons provide the dominant contributions. 

The total cross sections are calculated by integrating over the allowed kinematical 
region, accounting for the flnal state neutrino blocking. Figure ^ shows the energy 
dependence of the cross sections at the flducial values of hb = 0.64 fm~'^ and T = 20 
MeV. The neutrino chemical potentials in matter with and without hyperons are set by 
the equilibrium conditions. Neutrino flnal state blocking effectively suppresses the cross 
sections for E^,^ < fi^^. Neutrinos of higher energy are unaffected by blocking; hence, the 
cross section grows rapidly. In the degenerate regime, however, lepton number and energy 
diffusion are dominated by neutrinos close to the Fermi surface. (Initial state probabilities 
of high energy neutrinos are small.) Thus, an appropriate choice for the neutrino energy 



is the local neutrino chemical potential, i.e. E^^ = ^y^. In Figure |T0], we show the density 
dependence of the cross-sections (upper panel) in normal and strangeness-rich matter. The 
corresponding neutrino chemical potentials are shown in the lower panel. The presence 
of strangeness causes the neutrino concentrations to be signiflcantly larger than those in 
normal matter. As a result, the cross sections rise sharply with density in the presence of 
strangeness, chiefly due to the substantial increase in the neutrino energy. 



4.3. The Cooling Phase 

Here, we present results for thermal electron type neutrinos. The extension to /i and 
r type neutrinos and their anti-particles, by incorporating the appropriate neutrino-lepton 
coupling at tree level, is straightforward. (Neutrino coupling to leptons is flavour speciflc. 
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but their coupling at tree level to the baryons is not.) The transport of non-degenerate 
thermal /i and r neutrinos is important at all times, and thermal electron neutrino transport 
is important subsequent to deleptonization. 

In the cooling phase, neutrinos in the interior may be assumed to be in thermal 
equilibrium with Fermi-Dirac distributions and zero chemical potential; a typical thermal 
energy is E,j^ ~ vrT. Thus, a representative energy for thermal neutrinos lies in the range 



100-200 MeV. In Figure |TT], differential cross sections are shown at a baryon number density 
of ub = 0.64 fm^'^ for E^^ = 200 MeV. Scattering from the S~ hyperons dominates the 
contributions from the other particle species in matter at small u. 

The total cross-section per unit volume or the inverse collision mean free path is shown 



in Figure |I2| for E^^ = 200 MeV in matter with and without hyperons. The results highlight 
the density dependence of the cross sections. Here, since the electron neutrino chemical 
potential is zero, neutrino final state blocking is unimportant. The presence of hyperons 
in the cooling phase renders the interior more opaque to neutrinos relative to nucleons 
only matter. In fact, contributions from hyperons dominate the total cross section at high 
density. Clearly, composition can play an important role in neutrino scattering also during 
the cooling phase. 

Unlike in the deleptonization phase, where only neutrinos close to the Fermi surface 
contribute to diffusion of lepton number, the transport of energy in the cooling phase is 
determined by the cross sections of neutrinos of all energies. To understand the energy 
dependence, it is instructive to consider the contributions from negative uj (neutrino gains 



energy) and positive u (neutrino loses energy) separately (see Figure |l^). The cross section 
grows linearly for scattering when the neutrinos gain energy, while the cross section grows 
approximately as E^ for the case in which the neutrinos lose energy. It is important to 
note that, for thermal neutrinos, both contributions are important. The negative u phase 
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space dominates for low energy neutrino {E^, < ttT) scattering, while the positive cu phase 
space dominates at high energy {E^, > ttT) . The precise form of the energy dependence of 
the neutrino scattering rates is important in calculations of the energy averaged mean free 
paths for thermal neutrinos (Prakash et al., 1996). Our results suggest that an appropriate 
parametrization would be 



for the energy dependence for thermal neutrinos in the interior regions. 

A detailed discussion of the various averaging schemes employed in more complete 
treatments of neutrino transport is beyond the scope of this work. We hope to address 
these and related issues at a later time. 



Our aim here has been to elucidate the effects of composition and of strong interactions 
of the ambient matter on neutral current neutrino scattering cross sections during the 
deleptonization and cooling phases of the evolution of a newly born neutron star. Towards 
this end, we have calculated the neutrino scattering cross sections in protoneutron star 
matter, whose constituents exhibit varying degrees of degeneracy during the evolution of 
the star. In both phases, the composition of matter is chiefly determined by the nature 
of strong interactions and whether or not neutrinos are trapped in matter. In addition 
to the standard scenario, in which the strongly interacting particles are only nucleons, we 
have explored the influence of the possible presence of strangeness-bearing hyperons on 
the neutrino scattering cross sections. An important feature of our calculations is that the 
neutrino opacities are consistent with the equation of state of matter at finite temperature 
and density. 



V-^a{E) = A{nB, T)E + B{nB, T)E^ 



(40) 



5. 



CONCLUSIONS 
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We have identified neutral current neutrino interactions with hyperons that are 
important sources of opacity. Significant contribution to the neutrino opacity arises from 
scattering involving negatively charged E~ and S~ hyperons, chiefly due to their large 
vector couphngs, \Cv\ ~ 2. Although the contributions from the neutral A and E° are 
smaller than those of nucleons, these particles, when present, furnish baryon number which 
decreases the relative concentrations of nucleons. This leads to a larger opacity relative 
to nucleons only matter. The neutrino cross sections depend sensitively on the Fermi 
momenta and effective masses of the various particles present in matter. Whether or not 
a particular hyperon is present depends on the many-body description of charge neutral 
beta-equilibrated matter. We find that as long as one or the other hyperon is present, the 
cross sections are significantly modified from the case of nucleons only matter. 

In the deleptonization phase (lepton number fraction Yie = 0.4 and T ~ 20 — 30 
MeV), electron abundances are significantly larger than those in a cold catalyzed star, 
since neutrinos are trapped in matter, whether or not hyperons are present. Consequently, 
electrons contribute nearly as much as neutrons to the opacities. In neutrino trapped 
matter, the appearance of negatively charged hyperons (e.g., E~) is delayed to higher 
densities (relative to neutrino-free matter); also, their abundances are suppressed. However, 
the presence of neutral hyperons, such as the A, results in neutrino abundances that grow 
with density. This leads to significant enhancements in the cross sections for neutrinos (of 
characteristic energies close to the local neutrino chemical potential) compared to those in 
normal nucleonic matter. 

Modifications due to strangeness in the cooling phase are quantitatively different from 
those in the deleptonization phase. In the cooling phase, in which matter is nearly neutrino 
free, the response of the S~ hyperons to thermal neutrinos is the most significant. Although 
a comparable number of A hyperons are present, neutrinos couple weakly to this species. 
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and, hence, their contributions are significant only at high density. 

Our findings here suggest several directions for further study. The extension to include 
correlations between the different particles and RPA corrections to the results obtained 
here is under progress. The presence of charged particles, such as the E~, could make 
available low energy collective plasma modes through electromagnetic correlations, in 
addition to the scalar, vector and iso-vector correlations. Calculations of neutrino opacities 
from charged current reactions (which are important during the deleptonization phase), 
in strangeness-rich matter whose constituents exhibit varying degrees of degeneracy, are 
required for a complete description of the evolution. This will be taken up in a separate 
work. Effects of strangeness on lepton number and energy transport may be studied by 
employing energy averages (Rosseland means) of the opacities in present protoneutron 
codes. Useful tables of such average opacities will be made generally available. With new 
generation neutrino detectors capable of recording thousands of neutrino events, it may be 
possible to distinguish between different scenarios observationally. 
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APPENDIX A 

The integral form of the imaginary parts of the various polarizations that characterize 
the system's response to the neutrinos are collected here. The causal component of the 
density dependent polarizations in Eqs. (p!^-|T3|) are given by Saito et al (1989). For space 
like excitations (go < IqD and g^ < 0, they are given by 

Imn,(go,g-) = .A / ^ - \p\^ cos^ 6 

'^^ J (27r)3 E;E;^g ' ^ ^ 

r d'p ql/2 - \p\\l - cos' 



where 



Imnr(go,g) = -A / — 4^ ^ ^- e+, (2) 

2 J (27r)3 E;E;+g ' ^ ' 

Imn,(go,g-) = -vtA/^ ® ' 



e± = F±(i?;, [5(go - (i?;+, - i?;)) + siqo - (i?; - , (s) 
F±(E,E*) = /+(i5;)(i-/+(i5;*)±/_(i5;)(i-/_(E*), (6) 



e; = ^H2 + M*2. (7) 

(8) 

Above, A = 2 is the spin degeneracy factor. The particle distribution functions f±{E*) are 
the Fermi-Dirac distribution functions 

l + exp(^) 



where u is the effective chemical potential defined in Eq. (|30|). The angular integrals are 
performed by exploiting the delta functions, and the three dimensional integrals can be 
reduced to the following one dimensional integrals: 

A q] 



\ POO 

Imni(go,g-) = -^^j^ dE [{E + q,/2f - \q\' / A] 
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X [F+iE,E + qo) + F+iE + qo,E)], (10) 
TT |g| 

X [F+{E,E + qo) + F+{E + qo,E)], (11) 
A 



ImnT(go,g1 = ^fk\ dE [^E* + q,f + \q\^/A+\q\'M*'/ql)] 
87r a Je- 



Imn5(go,g) = --—AM*'-ql/A) dE [F+{E, E + qo) + F+{E + qo, E)] , (12) 
47r|g| Je_ 

AM* /""^ 

ImnM(go,g) = ^Tl dE [2E + qo][F-{E,E + qo) + F-{E + qo,E)]. (13) 

OTT O Je- 



The axial and the vector-axial polarizations entering the neutrino differential cross-sections 
are related to the scalar and mixed polarizations defined above through the following 
relations 

lmUAiqo,q) = -^^^j^—j^lm Usiqo, q) , (14) 

Im UvAiqo, q) = 2|g|2M* ^ • (^^) 

The causal polarizations are related to the retarded or time ordered polarizations through 

Im n^(go, q) = tanh (^^) Im n^(go, q) • (16) 

The phase space integrals 

In = tanh (^^) dE E- [F{E, E + q^) + F{E + go, E)] (17) 

simplify the representation of the retarded polarizations. The one dimensional phase 
space integrals are explicitly evaluated and expressed in terms of polylogaritmic functions. 
This representation is particularly useful in understanding the behavior with density and 
temperature and is given in Eq. (^Tj) in §2. For our purpose, only the three integrals Iq, Ii, 
and /s are required. Their analytical representations are also given in §2. 
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FIGURE CAPTIONS 

Fig. 1-Particle fractions, Fj = rii/nb at zero temperature, for the model GMl. The 
upper panel refers to nucleons-only matter. The lower panel shows results in matter with 
hyperons. 

Fig. 2-Particle fractions in the deleptonization phase (T = 20 MeV). The lepton fraction 
lie = + is chosen to be 0.4. 

Fig. 3-Particle fractions in the cooling phase (T = 50 MeV and Yy^ = 0). 

Fig. 4-Neutrino differential scattering cross sections in pure neutron matter at T = 10 McV. 
Curves labelled 1 and 2 refer to baryon densities of ub — 0.32 fm~^ and ub = 0.64 fm~^, 
respectively. Curves labelled 1(a) and 2(a) are results with interactions as in a mean field 
theoretical model. Those labelled 1(b) and 2(b) are results for a free Fermi gas. 

Fig. 5-Neutrino differential scattering cross sections in pure neutron matter at T = 10, 20, 
and 50 MeV. Interactions as in a mean field theoretical model are included. Upper panel 
results are for ub — 0.16 fm~^ and lower panel results are for ub — 0.64 fm~^. 

Fig. 6-Neutrino differential scattering cross sections in charge equilibrated matter at zero 
temperature. Upper panel results are for matter with nucleons only. The lower panel shows 
results in matter with hyperons. 
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Fig. 7-Same as in Figure 6, but for charge equilibrated matter at T = 30 MeV. 

Fig. 8-Neutrino differential scattering cross sections in charge equilibrated neutrino-trapped 
matter in the deleptonization phase. Upper panel results are for matter with nucleons only. 
The lower panel shows results in matter with hyperons. 

Fig. 9- Neutrino total scattering cross sections versus neutrino energy in charge equilibrated 
neutrino-trapped matter in the deleptonization phase. Blocking of final state neutrinos is 
included. 

Fig. 10-Upper panel: Neutrino total scattering cross sections in charge equihbrated 
neutrino-trapped matter in the deleptonization phase. The neutrino energy is set equal 
to the local neutrino chemical potential at each density. Lower panel: Electron neutrino 
chemical potential. 

Fig. 11-Neutrino differential scattering cross sections in charge equilibrated neutrino- free 
matter in the cooling phase. Upper panel results are for matter with nucleons only. The 
lower panel shows results in matter with hyperons. 

Fig. 12-Neutrino total scattering cross sections in charge equilibrated neutrino-free matter 
in the cooling phase. Important individual contributions to the total cross sections are also 
shown. 
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Fig. 13-Energy dependence of the neutrino total scattering cross sections in charge 
equihbrated neutrino-free matter in the coohng phase. Contributions from the positive 
energy transfer and the negative energy transfer processes, and their sum are shown in 
matter with nucleons only (upper panel) and in matter with hyperons (lower panel). 
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TABLE CAPTIONS 

TABLE 1. NOTE - Coupling constants derived assuming SU(3) symmetry and the 
constituent quark model for the hadrons. Numerical values are quoted using D=0.756 , 
F=0.477, sin2^vF=0.23 and sin^^ = 0.231 (Gaillard & Sauvage 1984). 

TABLE 2. NOTE.- Constants from Glendenning and Moszkowski (1991). 
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TABLE 1 

NEUTRAL CURRENT VECTOR AND AXIAL COUPLINGS 

Reaction Cy Ca 

Ue + e^Ue + e 1 + 4 sin^ Ow = 1.92 1 

i/e + /i ^ z^e + At -1 + 4 sin^ ew = -0.08 -1 

Ui + n^Ui + n -1 -D-F= -1.23 

Ui+p^Ui+p 1 -4sm2eH^ = 0.08 D + F = 1.23 

+ A^Ui + A -1 -F-D/3 = -0.73 

+ E-^i/j + E- -3 + 4sin2^H^ = -2.08 D - 3F = -0.68 

+ E+ ^ i/i + E+ 1-4 sin^ 9w = 0.08 D+F =1.23 

i/j + E° ^ i/i + E° -1 D-F^ 0.28 

i/i + 5- ^ i/i + 5- -3 + 4 sin^ % = -2.08 D - 3F = -0.68 

z/i + SO ^ i/i + SO -1 -D-F = -1.23 

i/i + E^^i/i + A 2i:)/y3 = 0.87 
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NUCLEON-MESON COUPLING CONSTANTS 



Model ^ 




9p_ 


b 


c 


M 


(fm) 


(fm) 


(fm) 









GMl 3.434 2.674 2.100 0.00295 -0.00107 0.70 
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Fig. 1.- 
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Fig. 3.- 
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Fig. 7.- 
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Fig. 9.— 
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Fig. 10.— 
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Fig. 11.— 
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Fig. 13.— 



